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Edge states in the three-quarter filled system, ai-(BEDT-TTF)2l3 
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We study the edge states in the two-dimensional conductor a- (BEDT-TTF) 2 13 theoretically. 
We show that the Dirac points and the edge states appear at the 3/4 and 1/4 filling as well 
as the half filling, due to four sites in the unit cell. This situation is in contract with the 
graphene, where the Dirac points and the edge states appear only at the half filling case. The 
edge states exist in both vertical and horizontal edges. For the 3/4 filled case it is shown that 
there exists edge states for all possible edges in the regions of \k y \ > \K y \ or \k y \ < \K y \ for the 
vertical edges, and \k x \ > \K X \ or \k x \ < \K X \ for the horizontal edges, where ±(K X , K y ) are 
the positions of the Dirac points in the bulk system, depending on the choice of edges. 
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1. Introduction 

The massless Dirac particles in the quasi-two- 
dimensional organic conductor cv-(BEDT-TTF) 2 l3 un- 
der pressure 1 - 1 have been predicted theoretically by 
Katayama, Kobayashi and Suzumura, 2 -* and it has been 
confirmed by the interlayer magnetoresistance. 3,4 ) 

The energy dispersion of the massless Dirac particles 

En±v\k-k \, (1) 

has been realized in graphene. 5, 6 ) Two bands touch at 
the Dirac points ±kn. The edge states exist in graphene 
and have been studied by many authors. 7 ~ 10 ) 

Graphene has a honeycomb lattice structure and it has 
two sites in the unit cell. There exists one electron per 
each site, i.e. the band is half-filled in graphene. On the 
other hand, a-(BEDT-TTF) 2 I 3 has four sites in the unit 
cell as shown in Fig. 1, where the small transfers between 
planes are neglected. The band made from the molecular 




Fig. 1. (Color online) The schematic figure of two-dimensional 
lattice for a-(BEDT-TTF) 2 I 3 . There are four sites (1, 2, 3,and 
4) in the unit cell, three kinds of intra chain transfers (ai, 02 
and 03), and four kinds of inter-chain transfer integrals (bi, &2, 
63 and 64). 



orbits of BEDT-TTF molecules is 3/4 filled, since one 
electron is moved from two BEDT-TTF molecules to I3 
molecule. The four sites in the unit cell make the Dirac 
points appear at 1/4 and 3/4 filled cases besides the half- 
filled case. 

We study the 3/4 filled band in a-(BEDT-TTF) 2 I 3 
by using the tight-binding model in two dimension. The 
bulk properties of ev-(BEDT-TTF) 2 l3 can be expressed 
by the effective two-band model, n ~ 13 ) which is sim- 
ilar to the model of the graphene with next-nearest 
hoppings. 12, 14 ' 15 -* However, when we consider the edge 
states, the effective two-band model cannot be applied. 

In this paper, we study the edge states in a-(BEDT- 
TTF) 2 i3 at 3/4 filling and compare these edge states 
with the edge states in graphene at half-filling. We show 
that the edge states always appear either in the region 
of \k y \ > \K y \ or \k y \ < \K y \ for the vertical edges and in 
the region of | A;^ | > \K X \ or \k x \ < \K X \ for the horizon- 
tal edges, depending on the choice of the possible edges, 
where k x are k y are the wave number of the eigenstate 
and ±(K X , K y ) are the wave number of the two Dirac 
points in the bulk system. 

2. Model 

We study the tight-binding model for the quasi-two- 
dimensional conductor a-(BEDT-TTF) 2 I 3 . There are 
four sites in the unit cell, which we label as sites 1, 2, 

3, and 4, as shown in Fig. 1. There are the hoppings 
along the y direction, a%, a 2 and 0,3, and the inter-chain 
hoppings b\ — 64. The system has the inversion symme- 
£ r y_i6-i8) rp ne mvers j on center is located at the site 3, 
the site 4, and the center of the sites 1 and 2. With the 
inversion the sites 1 and 2 are exchanged each other but 
the sites 3 and 4 remain unchanged. The sites 1, 2, 3 and 
4 are also called AI, All, B and C, or A, A', B, and C, 
respectively. We also take account the site energies ei, 
e 2 , £3 and €4. When t\ = e 2 , the inversion symmetry is 
conserved. 18 - 1 The energy of electrons in the tight binding 
approximation is given by the equation 
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Fig. 2. (Color online) The Dirac cone of a-(BEDT-TTF) 2 I 3 . 
The upper two bands touch at the Dirac points, ±(K X , K y ) = 
±(0.398tt, -0.353tt) 
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where ipn/ T 



ipr^Jm are the wave functions at sites 1 — 
4 in the unit cell (n, m) with integers n and m. The 
parameters are estimated by Kondo et al. 18,19 - 1 and used 
by Kobayashi et al. 20) as ax = -0.028 - 0.0025p, a 2 = 
0.048 + 0.008p, a 3 = -0.020 + 0.0005p, 6i = 0.123, 6 2 = 
0.140 4- 0.0015p, 6 3 = -0.062 - 0.002p and 6 4 = -0.025, 
where p is the uniaxial strain in the y direction. In this 
paper we use the parameter at p — 4 kbar, i.e. a\ — 
-0.038, a 2 = 0.08, a 3 = -0.018, b x = 0.123, b 2 = 0.146, 
6 3 = -0.07 and 6 4 = -0.025. 



3. Bulk system 

If the system is infinite or periodic with respect to both 
x and y directions, the eigenstates are written as 
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Then the energy is obtained by 
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If ax = a 2 = a 3 = (i.e. Cx = C 2 = 0) and £1=62 = 
63 = £4 = 0, the eigenvalues of the matrix in Eq. (9) are 
obtained by Mori 17 ) as 



E k = ±\ \Dx\i + \D 2 \i ± J(D{ + Di)(Dl 



D* 



(14) 

The condition for the Dirac points at the 1/4 and 3/4 
filled band is obtained for this case of Cx = C 2 = as, 



Dx = ±iD 2 . 

This condition is fulfilled at (k x ,k y ) — 
K x and K y is given by 

b 2 e lK * + b 3 



(15) 

(K x , K y ), where 



bxe m + 04 

In order to have real solutions, the absolute value of the 
right hand side in Eq. (16) should be unity, from which 
we obtain 

b\ 



COS if T 



b\ 



b\ 



b\ 



(17) 



2(6263-6164) ' 

which has been obtained by Mori. 17 ) 

When Cx and C 2 are not zero, which is the case in 
a-(BEDT-TTF) 2 l3, the eigenvalues are complicated, al- 
though the analytical expression is possible as a solution 
of the quartic equation. 

In Fig. 2, we show the 3D plot of the energy with 
the parameters for a-(BEDT-TTF) 2 I 3 with p = 4 kbar 
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and £1 = 62 = £3 = £4 = 0. In this figure we show 
the fourth and the third bands from the bottom of the 
energy (i.e. the first and the second band from the top 
of the energy) and the part of the second band from the 
bottom of the energy (the third band from the top of 
the energy). The first and the second bands from the top 
of the energy touch at two Dirac points ±(K X , K y ) = 
±(0.3987r, — 0.3537t). In this choice of parameters a finite 
gap exists between the second and the third bands as 
shown in Fig. 2. In this paper, however, we focus on the 
edge states in the three-quarter filled band and do not 
study the half-filled case, since we are interested in a- 
(BEDT-TTF) 2 I 3 , which is the system with the 3/4-filled 
band. 

4. Edge states 

4-1 vertical edges 

In this section we study the system with edges. The 
canted edges may be possible, but we consider the verti- 
cal and horizontal edges in this paper. 

First we consider the vertical edges. We study two pos- 
sibilities for each edge. The left edge and the right edge 
can be either the chain of the sites 1 and 2, or the chain 
of the sites 3 and 4 (see Fig. 1). Consider the ribbon of 
2 x L chains (where L is integer) in which the left edge 
is the chain of sites 1 and 2, and the right edge are the 
chain of the sites 3 and 4. We call this system as the (12- 
34) edge. The other choices of the edges we will study 
for the vertical edges are the boundaries with the chain 
of the sites 3 and 4 at the left edge and the chain of the 
sites 1 and 2 at the right edge, which we call the (34-12) 
edge, the boundaries with chains of the sites 1 and 2 at 
both edges (the (12-12) edge), and the boundaries with 
chains of the sites 3 and 4 and at both edges, (the (34-34) 
edge). Note that the (34-12) edge has 2L chains, while 
the (12-12) edge and the (34-34) edge have 2L + 1 chains. 

We assume the periodic boundary conditions in the 
y direction for the systems with vertical edges, which 
means that we study the very long vertical ribbon or the 
tube similar to carbon nanotube. 

Similar to the honeycomb lattice, 10 ) we perform the 
Fourier transformation with respect to y as 
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The energy and the eigenstate are obtained by the equa- 
tion 
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where 0, M, N, are 4x4 matrix given by 
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and * is the vector with 4L wave functions, 



(i) 

ky 

(2) 

(3)" 
n,k y 



L,k y 
L,k y 



(29) 



L.k v 



J 



In Eq. (19), we had introduced the parameter a, which 
describe the boundary conditions. For the open boundary 
condition, in which the edge states may exist, we should 
take a = 0. The periodic boundary conditions can be 
obtained by taking a = 1. 

In the same way we take the boundary with the 3 and 
4 sites at the left edge and 1 and 2 sites at the right 
edge, which we call (34-12) edge. The energy of the (34- 
12) edge is the same as that of the (12-34) edge except 
that the left and the right edge states are exchanged. 

In the (12-12) edge and the (34-34) edge, we cannot 
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apply the periodic boundary conditions in the x direc- 
tion. The matrix size is (4L + 2) x (4L + 2) in these edges. 

Note that the (12-12) edge and the (34-34) edges are 
symmetric with respect to inversion, while the (12-34) 
edge and the (34-12) edge are not symmetric but they 
are exchanged each other by inversion. 

4-2 horizontal edges 

For the horizontal edges, the lower and upper bound- 
aries have the zigzag shape. Each edge consists with one 
pair of the four possibilities, i.e., sites 4 and 2, 2 and 3, 3 
and 1, or 1 and 4. There are 16 possibilities for the choice 
of the horizontal edges. We call these 16 possible edges 
as the (42-42) edge, the (42-23) edge, and so on. Consid- 
ering that the sites 1 and 2 are exchanged each other by 
inversion while sites 3 and 4 are not changed by inver- 
sion, we obtain that the (42-14) edge, the (23-31) edge, 
the (31-23) edge and the (14-42) edge are symmetric with 
respect to inversion. The (31-42) edge, for example, be- 
comes the (14-23) edge by inversion. For the horizontal 
edges, we can perform the Fourier transformation with 
respect to x and the energy is labeled by k x . 

There are 4L eigenstates for each k x in the (42-31) 
edge, the (23-14) edge, the (31-42) edge and the (14-23) 
edge, where integer L is the number of each sites (1—4) in 
each vertical chain. In these four cases we can apply the 
periodic boundary conditions with respect to y direction, 
as in the cases of the (12-34) edge and the (34-12) edge 
for the vertical edges. If the periodic boundary conditions 
with respect to k y direction are taken, the above four 
cases give the same energy dispersion. 

There are 4L + 1 eigenstates for the (42-14) edge, the 
(23-42) edge, (31-23) edge and the (14-31) edge, where 
the lowermost and the uppermost sites are the same. 
When the lowermost and the second lower sites are the 
same as the second upper and the uppermost sites, such 
as the (42-42) edge, the (23-23) edge, the (31-31) edge 
and the (14-14) edge, there are 4L + 2 eigenstates for 
each k x . When the second lower sites are the same as the 
second upper sites, such as the (42-23) edge, the (23-31) 
edge, the (31-14) edge and the (14-42) edge, there are 
4L + 3 eigenstates for each k x . 

5. Results 

There arc N = 4 x L eigenstates and energies for each 
k y for the (12-34) edge and the (34-12) edge. We plot 
the energy as a function of k y in Fig. 3, where we use 
different colors for each band of L states. There are edge 
states at the 1/4, 1/2 and 3/4 fillings of the band. In this 
paper we focus on the edge state at 3/4 filling case, since 
3/4 filling is realized in a-(BEDT-TTF) 2 I 3 . 

In Fig. 4 we plot the energies near the 3/4 filling for 
the system with periodic boundary conditions for the 
x-direction, the (12-34) edge, the (12-12) edge and the 
(34-34) edge. If the system is periodic with respect to x 
i.e., a — 1, we obtain the projection of the 3D plot of the 
energy (Fig. 2) as shown in Fig. 4(a). 

The eigenstates of the (12-34) edge at k y = 0.6tt and 
k y = for the L-th state from the top and that at k y = 
0.67T for the (L + l)-th state are shown in Fig. 5 (a) , (b) 



(12 -34) edge L=25, N=4 L 
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4 : b = 0.123, b 2 = 0.146, t> 3 = - 0.07, b = - 0.025 
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Fig. 3. (Color online) The energy spectrum in the system with 
four sites in the unit cell with the (12 - 34) edge. We take the 
parameters for the a-(BEDT-TTF)2l3 under the uniaxial pres- 
sure 4 kbar in the a-direction. There are four bands, each of 
which contains the L eigenvalues. The different brightness of lines 
(purple, green, red, and blue in the online version) represents the 
different bands. Thick lines with left triangles are the edge states 
on the left edge, and thick lines with right triangles are the edge 
states on the right edge. 

and (c), respectively. The L-th state from the top, i.e. the 
bottom of the highest band, is the edge state, which is 
localized at the left edge for \k y \ = 0.67T > \K y \ « 0.3537T 
(Fig. 5 (a)) and at the right edge for \k y \ = < \K y \ 
(Fig. 5 (b)). We define the localization length, no, as 

l< ) fcB | 2 «exp(-n/n ), (30) 

for the edge state localized at the left edge and 

\^ k f <x exp(-(L - n)/no), (31) 

for the edge state localized at the right edge. It is ob- 
tained that no « 1.9 at k y — 0.67T and n$ w 0.37 at 
k y = 0, as seen in Fig. 5(a) and (b). Other states such as 
the (L + l)-th state, for example, are not the edge state 
as seen in Fig. 5(c). 

The edge states exist only at \k y \ > \K y \ for the (12- 
12) edge and at \k y \ < \K y \ for the (34-34) edge, as seen 
in Fig. 4 (c) and (d). In these cases the edge states are lo- 
calized in both left and right edges. In general, when the 
localization length, no, of the edge state is much smaller 
than the width L (n <C L), the edge states at each edge 
can be treated as independent states, and these states de- 
generate due to inversion symmetry. Otherwise, the edge 
states at each edge interact each other and the degen- 
eracy is lifted. Then the "bonding" and "anti-bonding" 
states of the each edge states become the eigenstates. In 
our numerical calculation for L = 100, which is much 
larger than the localization length, the L-th and the 
(L + l)-th states from the top have the same energy 
within the numerical accuracy and the eigenstates are 
any linear combinations of the left and the right eigen- 
states, as shown in Fig. 6, where L-th and (L + l)-th 
states are localized at both edges. 
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Fig. 5. (Color online) Edge states near 2 filling for systems with 
(12 — 34) edges with L = 100. The squares of the absolute value 
at (a) k y = 0.67T and (b) k y = of the L-th state from 

the top, and at (c) k y = 0.6ir of (L + l)-th state from the top are 

-,(2) 



plotted as a function of n. The values of \*&„ a^l 2 : 6irl 2 



and l^nfi.eJ 2 in ( a ) and the valueS of l*n,0.6^| a and l*n,0.6,rl 
in (c) are almost the same in this scale of the figure. The expo- 
nentially localized character is seen in (a) and (b). 



A2) 



Fig. 4. (Color online) Energy spectrum near | filling for systems 
with periodic boundary and with edges. The parameters are the 
same as in Fig. 3. 
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Fig. 6. (Color online) Edge states near ^ filling for systems with 
(12 — 12) edges with L = 100. The squares of the absolute value 

*n fc„ = °- 67r ) of the L ~ th ( a ) and ( L + 1 )~ th states from the 

top are plotted as a function of n. The values of V&^q Gn ^^fo g^- 

and ^^ 4 q 67r are almost the same in this scale of the figure. The 
exponentially localized character is seen. 

From the above results we conclude that if the left or 
right edge is the chain with sites 1 and 2, the edge states 
exists in the region \k y \ > \K y \ near 3/4 filling. If the left 
or right edge is the chain with sites 3 and 4, the edge 
states exists in the region \k y \ < \K y \ near 3/4 filling. 
We summarize the existence of the edge states in Table I 
(a). 
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Fig. 7. (Color online) Energy spectrum of a-(BEDT-TTF) 2 l3 
with periodic boundary conditions as a function of the wave 
number k x . The parameters are the same as in Fig. 3. 

Next we study the horizontal edges. In Fig. 7 we plot 
the energy of the system with periodic boundary condi- 
tions as a function of k x . For various types of horizontal 
edges, we find that edge states always exist. We plot some 
examples of energy of the horizontal edges as a function 
of k x in Fig. 8. 

We obtain that in the 3/4-filled case there are the edge 
states localized near the lower edge of the (14-xy) and 
(31-xy) edges when \k x \ > \K X \, where xy is either 14, 
31, 23, or 42. Since (14-xy) and (31-xy) edges are changed 
to (y'x'-42) and (y'x'-23) edges by inversion, where y'x' 
is either 42, 23, 31, or 14, respectively, (note that sites 1 
and 2 are exchanged each other by inversion, while sites 
3 and 4 are not changed.), there are edge states localized 
near the upper edge of (y'x'-42) and (y'x'-23) edges when 
|*x| > \K X \. 

We also obtain that there are the edge states local- 
ized near the lower edge of the (23-xy) and (42-xy) edges 
when \k x \ < \K X \, and there are the edge states localized 
near the upper edge of the (y'x'-31) and (y'x'-14) edges 
when \k x \ < \K X \. Table I (b) shows the summary of the 
existing region of the edge states for the horizontal edges. 

The edge states in the same rows in table I (b), for ex- 
ample, the lower 14 and upper 42, have the same energy 
dispersion. The edge states in different rows in table I (b) , 



(a) vertical edge 



left edge state 


right 


edge state 


12 \ky\ > \Ky\ 


12 


\ky\ > \Ky\ 


34 \k y \ < \Ky\ 


34 


\ky\ < \Ky\ 


(b) horizontal edge 


lower edge state 


upper 


edge state 


14 \k x \ > \K X \ 


42 




31 1^1 > \K X \ 


23 


&x I^x 


23 \k x \ < \K X \ 


31 


| fcx | *C | K x | 


42 Ifc^ < \K X \ 


14 


I &x < K x | 



Table I. Existing region of the edge states in «-(BEDT-TTF) 2 I 3 
at 3/4 filling. See Fig. 11(a). 



vertical edge 



left, right 


edge state 


zigzag 

bearded (Klein's edge 21 ') 


\ky\ > \Ky\ 
\ky \ < \Ky\ 


horizontal edge 


lower, upper 


edge state 


armchair with Klein's edge 22 ' 
armchair 


\k x \ > \K X \ 
\k x \ < \K X \ 



Table II. Existing region of the edge states in graphene. For the 
isotropic honeycomb lattice, \K X \ = and \K y \ = ^tt. Then no 
edge states exist in the armchair edge. However, edge states in 
armchair edge exist in the anisotropic case, 10 ' where \K X \ 0. 
See Fig. 11(b) and (c). 
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(a) 



(31-42) edge, (14-23) edge L=30, N=4 L 
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(42-42) edge, (14-14) edge L =3(), N=4 L+2 




L=30, N=4 L+3 




(42-31) edge, (23-14) edge L=30, N=4 L 




Fig. 8. (Color online) Energies near the 3/4 filling in some hori- 
zontal edges. The energies of the edge state at \K\ > \K X \ near 
the lower edge of sites 1 and 4 or near the upper edge of sites 4 
and 2 are close to the energies of the top of the second band, as 
shown in (a), (b), (c) and the inset of (c), which is the magnified 
figure. 
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Fig. 9. (Color online) Wave function of the (a) (L + l)-th, (b) 
(L + 2)-th, and (c) (L + 3)-th states from the top of the energy 
for the (14-42) edge at k y = 0.8n. We take the same parameters 
parameters as in Fig. 3. The (Z/ + l)-th and (L + 2)-th states are 
edge states, but (L + 3)-th state is not the edge state. 



such as the lower 31 and upper 42, have different energies 
as a function of k x . In Fig. 8(a), we plot the energy as a 
function of the k x in the (31-42) edge (the (14-23) edge 
has the same energy). The L-th and the (Z + l)-th states 
from the top of the energy are edge states at \k x \ > \K X \. 
The edge states in the (L + l)-th state (at \k x \ > \K X \) 
has a relatively large localization length and the energy 
of the edge state is close to the second band. The energy 
of the [L + l)-th state at \k x \ > \K X \ in the (31-42) edge 
(Fig. 8 (a)) is same as the (L + l)-th state at \k x \ > \K X \ 



8 J. Phys. Soc. Jpn. 



Full Paper 



Y. Hascgawa and K. Kishigi 



in the (42-42) edge (Fig. 8 (b)) and the (L + l)-the and 
the (L + 2)-th states in the (14-42) edge (Fig. 8 (c) and 
the inset). In Fig. 9 we plot the wave function of the 
(L + l)-thc, the (L + 2)-th and (L + 3)-th states in the 
(14-42) edge at k x — 0.87T as a function of the site num- 
ber to. It is clearly seen that the (L + l)-th state and 
the (L + 2)-th state are the edge states localized at both 
lower and upper edges, while the (L + 2)-th state is not 
the edge state. 

The energies as a function of k x for the edge states of 
the (23-xy) edge and the (42-xy) edge are different, as 
seen two curves between the first and the second bands 
(L-ih and (L+l)-th states) from the top of the energy at 
\k x \ < \K X \ in Fig. 8 (d). As summarized in Table I, all 
edges have the edge states at \k y \ > \K y \ or \k y \ < \K y \ 
for the vertical edge and at | fc^ | > \K X \ or | A;^ | < \K X \ 
for the horizontal edge. 

The edge states in a-(BEDT-TTF) 2 I 3 obtained above 
are compared to the edge states in graphene with 
anisotropic hoppings. In graphene there are two sites in 
the unit cell, forming A and B sublattices. These sites are 
exchanged by inversion. The edges in the graphene are 
classified as zigzag, bearded, zigzag-bearded and arm- 
chair. 10 ) We can also consider the armchair with the 
Klein's edge. 22 ) These edges are shown in Fig. 10. 

We consider the anisotropic cases with different hop- 
pings between nearest neighbors in different direc- 
tions. 10 - ) For the isotropic honeycomb lattice the edge 
states exist in the zigzag, bearded and zigzag-bearded 
edges for \k y \ > 27r/3, \k y \ < 27r/3 and any k yi respec- 
tively, while there is no edge states in the armchair edges. 
The Dirac points in the isotropic honeycomb lattice are 
located at (0, ±47r/(3a)), where a is a lattice constant of 
the honeycomb lattice, as shown in Fig. 11(b). If we con- 
sider the anisotropic cases, however, it has been shown 
that the edge states exit even in the armchair edges. 10 ) 
The reason for the absence of the edge states in the 
isotropic honeycomb lattice with armchair edges can be 
understood as follows. In the armchair edges, the edge 
states can exist in the region \k x \ < \K X \ where ±K X are 
the x component of two Dirac points. In the case of the 
isotropic honeycomb lattice, the projections of the two 
Dirac points into k x axis coincide at K x = 0. In that case 
there are no region of k x , which satisfies k x < \K X \ = 0, 
and as a result there are no edge states in the isotropic 
armchair edges. On the other hand, if we consider the 
anisotropic honeycomb lattice with K x ^ 0, we have the 
edge states at \k x \ < \K X \. Recently, it is also shown that 
the edge states exist in the isotropic honeycomb lattice 
with armchair edges moderated by the Klein's defects. 22 ) 
In Table II, the exiting region of the edge states are given. 

We speculate that if we can make the tilted edges with 
the tilting angle 9 = — &rct'A\\(K x / K y ), then coincidence 
of the projections of two Dirac points in the wave number 
along the tilted direction, as shown in Fig. 11(a). In that 
case we obtain the similar situation as in the isotropic 
honeycomb lattice with armchair edges and we expect 
that the edge states cannot exist. 

One of the differences between the edge states in a- 
(BEDT-TTF) 2 l3 and the edge states in graphene is that 



armchair with Kleins' s edge 




armchair 



Fig. 10. (Color online) Possible edges for the graphene (honey- 
comb lattice), zigzag, bearded, armchair and armchair with the 
Klein's edge. 

all four components of the wave functions decays expo- 
nentially away from the edge in edge states in a-(BEDT- 
TTF) 2 l3, while one component of the wave function is 
always zero and the other component decays exponen- 
tially in the edge states in graphene. The other differ- 
ence is the wave number dependence of the energy of the 
edge state. The energy of the edge states in a-(BEDT- 
TTF) 2 i3 has strong k y or k x dependence, which is in 
contrast with the edge states in graphene, where the en- 
ergy of the edge states is zero if the next-nearest-neighbor 
hoppings are not taken into account. The edge states in 
graphene have the k y dependent energy only when the 
next-nearest-neighbor hoppings are finite. 8 ' 9 ' 

The energy dispersion of the edge states results in the 
occupation of the edge states at \k y \ < \K y \ in (12-34) 
edge when the system is 3/4 filled, since the energy of the 
edge states is lower than the energy at the Dirac points, 
as seen in Fig. 4 (b). In that case only the edge states at 
the right edge are occupied. Similarly, the energy of the 
edge states for the horizontal edge have a wave-number 
dependence, as seen in Fig. 8. 

6. Site energy 

In this section we study the effect of the site energies 
ei. Note that 63 and 64 do not violate the inversion sym- 
metry but ei and e 2 violate the inversion symmetry if 
€\ 7^ e 2 . When ei = e 2 = e 3 = £4, the dispersion of the 
eigenstates is only shifted by the constant energy. There- 
fore, if we consider the system with inversion symmetry, 
£3 and €4 are the relevant parameters. Either ei or e 2 can 
be taken as a parameter for the breaking of the inversion 
symmetry. 

The energy gap is opened when inversion symmetry is 
broken by e\ or e 2 , as shown in Fig. 12. Even if the gap 
is opened by e\ or e 2 , the edge states still exist. As seen 
in Fig. 12 (a) the L-th state from the top in the system 
with (12-34) edge is the left edge state for \k y \ > K y and 
the right edge state for \k y \ < K yi which is the same as in 
the system with inversion symmetry (see Fig. 4(b)). The 
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Fig. 11. (Color online) The Dirac points for electrons on (a) 
o-(BEDT-TTF) 2 l3, (b) isotropic honeycomb lattice, and (c) 
anisotropic honeycomb lattice. The small striped circles are T 
points. The Dirac points are shown by blue open circles and red 
filled circles. Hexagons in (b) and (c) are the Brillouin zone. The 
green rectangles in (b) and (c) have the same areas as hexagons. 
If the hoppings between the nearest sites are anisotropic, the 
Dirac points move from the corners of the Brillouin zone as shown 
in (b) and (c). If the tilted edge in a- (BEDT-TTF) 2 13 can be 
made, the projections of the Dirac points into the direction of 



(12-34) edge 



L=100, N=4L 



0.21 




(b) 



(12-12) edge 



L=100, N=4L+2 




Fig. 12. (Color online) Energy spectrum near j filling for systems 
with (12-34) and (12-12) edges. The parameters are the same 
as in Fig. 3 except for ei = 0.02, which violate the inversion 
symmetry. 



(12-12) edge has the edge states for \k y \ > K y as shown 
in Fig. 12 (b). Because the inversion symmetry is broken, 
the left and right edge states are not degenerated. The 
L-th state from the top of the energy is the edge state 
localized at the left edge, and the (L + l)-th state is the 
edge state at the right edge. The localization length of 
these edge states are different as shown in Fig. 13 (a) 
and (b). 

As seen in Fig. 12, the edge states with t\ ^ have 
the similar strong k y dependent energy as in the system 
with inversion symmetry. As a result, if the system is 
3/4 filled, i.e. the chemical potential is in the middle of 
the gap of the bulk system (for example fi — 0.175 in 
Fig. 12 (a) and (b)), the localized state at the right edge 
is partially filled. 

The first and the second bands from the top of the 
energy still touch each other at the Dirac points even 
when we take finite 63 and £4 (as an example we take 
£3 = 0.1 and £4 = —0.1), as shown in Fig. 14 (a) for 
the periodic boundary conditions. In that parameters 
(£3 = 0.1 and £4 = —0.1), the edge states exist as in 
the case of £3 = £4 = (see Fig. 4 and Fig. 14). The 
localization length, however, differs from that in the case 
of £3 = £4 = 0. While the L-th state from the top at 
k y = in the (12-34) edge is the edge state localized 
at the right edge with the localization length n w 0.83 
(Fig. 15), the L-th state at k y = 0.87T is the edge state at 
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Fig. 13. (Color online) Wave functions of the L-th and (L + l)-th 
states from the top at k y = 0.8tt for systems with (12-12) edge. 
The parameters are the same as in Fig. 3 except for ei = 0.02, 
which violate the inversion symmetry. The L-th state from the 
top is localized at the left edge with the localization length no ~ 
2.90 and the (L + l)-th state from the top is localized at the right 
edge with the localization length no ~ 1.80. 

the left edge with the large localization length uq w 20 
(Fig 16). The large localization length is more clearly- 
seen in the (12-12) edge. The L-th and (L + l)-th states 
at k y — 0.87T are almost degenerate in energy in the sys- 
tem with L = 100. However the wave function of the 
L-th state are different from that of the (L + l)-th state 
as shown in Fig. 17. The square of the absolute values 
of the wave functions, fc | 2 , of L-th and (L + l)-th 
states from the top at k y = 0.8tt are asymmetric and 
symmetric with respect to n and they look like "anti- 
bonding" and "bonding" states, respectively. This means 
that these states are indeed the localized edge states but 
the system size (L = 100) is not large enough for the edge 
states to be treated as independent states at each edges. 
In these cases the estimation of the localization length 
by Eqs. (30) or (31) may have an ambiguity and should 
be taken carefully. The localization lengths for fc | 2 
should be same for all i but depend on the component 
i as we give the estimation in the figure captions. Even 
in these cases, the existing regions of the edge states are 
given in Table I with K y varied according to the choice 
of e3 and £4. 
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Fig. 14. (Color online) Energy spectrum near ^ filling for systems 
with periodic boundary and with edges. The parameters are £3 = 
0.1 and e 4 = -0.1. 
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Fig. 15. (Color online) Edge states near ^ filling for systems with 
(12 — 34) edges with L = 100. The squares of the absolute value 
fe = 0) °f ^ ne ^~th state from the top is plotted as a 

function of n. It is seen that this state is an edge state localized 
at the right edge with small localization length no ~ 0.83. 
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Fig. 16. (Color online) Edge states near ^ filling for systems with 
(12 — 34) edges with L = 100. The squares of the absolute value 
fe = 0) °^ ^ e -^"t n state from the top is plotted as a 

function of n. It is seen that this state is an edge state localized 
at the left edge with the localization length no m 20.4. 20.3, 
21.5, and 20.3 for i. | with i = 1 — 4, respectively,, which are 
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Fig. 17. (Color online) Edge states near ^ filling for systems with 
(12 — 12) edges with L = 100. The squares of the absolute value 
*W fc (k y = 0.6?r) of the L-th (a) and (L + l)-th (b) states 
from the top are plotted as a function of n. The exponentially 
localized character is seen. The localization lengths are estimated 
as n 17.4, 15.8, 50.0 and 16.4 for of the L-th state from 

the top with i = 1 — 4, respectively, and n « 22.8, 24.4, 14.2 
and 23.7 for these of the (L + l)-th state. 



obtained from the least square fit of 1 



for 5 < n < 20. 
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7. Summary and discussions 

In this paper we study the edge states in a-(BEDT- 
TTF)2l3 in the absence of the magnetic field, theoreti- 
cally. We have shown that the edge states exist in the 3/4 
filled band in the system with four sites in the unit cell, 
such as a- (BEDT-TTF) 2I3 . We study the vertical and 
horizontal edges. We show that all edges have the edge 
states as summarized in Table I. In the edge states all 
components of the wave function decays exponentially 
as a function of the distance from the edge. Since the 
edge state has the wave-number dependent energy, it is 
possible to occupy only the left edge state or the lower 
edge states at 3/4 filling system. For example, in the (12- 
34) or (34-12) edges only the localized state at the right 
edge or the left edge is occupied in 3/4 filling system, 
respectively. 

We also study the effect of the site energies, e\, £3 
and €4. When the inversion symmetry is broken by t\ ^ 
or e 2 7^ 0, a finite gap appears at the Dirac points. Even 
in that case the edge states exist as in the system with 
inversion symmetry. Due to the k y dependence of the 
edge states, the energy of the edge states can locate in the 
middle of the bulk band gap. Therefore, the insulator at 
the bulk sample has the partially filled edge states. These 
states are not the topologically protected edge states, 
which are recently predicted 23-25 ) and observed 26-28 ) in 
the topological insulators. 

The opening of the gap at the Dirac points can be 
possible either by the breaking of the inversion symme- 
try (in the honeycomb lattice with next-nearest-neighbor 
hoppings the necessary symmetry for the zero gap is 
not the inversion symmetry but the "averaged inver- 
sion symmetry" . 14 ' 15 )) or by the merging of the Dirac 
points. 11,29,30 ^ If the inversion symmetry is not bro- 
ken in (BEDT-TTF) 2 I 3 , the opening of the gap at the 
Dirac points occurs only when two Dirac points merge at 
k = k modulo a reciprocal vector, i.e. (k x , k y ) — (0, 0), 
(0, 7r), (71", 0) or (71", 7r), which is predicted to occur at high 
pressure. 11 ^ In the present choice of parameters, the fi- 
nite gap at the half-filling as seen in Fig. 2 is caused 
by the merging of the Dirac points at (7r,0). The edge 
states also exist even when the gap is opened either by 
the breaking of the inversion symmetry or the merging 
of the Dirac points. 

With finite £3 and £4 the gap remains zero at the Dirac 
points, but the localization lengths of the edge states may 
be changed drastically. The edge states with the large 
localization length can be realized if we can tune the site 
energies. 

The experimental observation of the edge states in a- 
(BEDT-TTF) 2 l3 will be possible and provide us more 



insight about the massless Dirac particles realized in the 
quasi- two-dimensional conductors. 
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